SU(4) composite fermions in graphene: New fractional quantum Hall states 
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Theoretical studies of the fractional quantum Hall effect (FQHE) in graphene have so far focused 
on the plausibility and stability of the previously known FQHE states for the interaction matrix 
elements appropriate for graphene. We consider FQHE for SU(4) symmetry, as appropriate for the 
situation when all four spin and valley Landau bands are degenerate, and predict new FQHE states 
that have no analog in GaAs. These result from an essential interplay between the two-fold spin 
and valley degeneracies at fractions of the form v — n/(2pn± 1), for n > 3. Conditions are outlined 
for the observation of these states and quantum phase transitions between them; the structure of 
these states and their excitations is also described. 
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I. INTRODUCTION 

Soon after the first experimental realization of 
graphene, a single-layer hexagonal form of carbon, the 
integral quantum Hall effect (IQHE)i was observed by 
Novoselov et al. and Zhang et ai£ at filling factors 
v n = 4n — 2. The Ai/ — 4 period is understood straight- 
forwardly as a consequence of the four-fold (near) de- 
generacy of Landau levels (LLs) in graphene; the Zee- 
man energy is small in comparison to the interaction en- 
ergy scale, and the pseudospin degree of freedom, which 
represents the two inequivalent Dirac-cones at the cor- 
ners of the Brillouin zone, does not couple to external 
fields if the two sublattices are equivalent. The offset 
—2 is a consequence of a Dirac-like effective Hamilto- 
nian, which produces a linear dispersion for low energy 
electronic states^. Recently, Zhang et alA have observed 
quantum Hall plateaus at v = ±4, ±1, 0. Although exter- 
nal effects can surely lift the four-fold degeneracy of the 
low-lying LLs, Yang, Das Sarma and MacDonald 5 have 
shown that the exchange interaction can break the SU(4) 
symmetry spontaneously, and the charged excitations are 
skyrmions& in the \n\ < 3 Landau levels at v = 1 (or, by 
particle-hole symmetry, at v — 3) where v = v— v n is the 
filling factor within the Landau level in question. This 
has been confirmed^ by exact diagonalization for n < 2 
with N < 9 particles. 

Although the fractional quantum Hall effect 8 (FQHE) 
has not yet been observed in graphene, it has been ex- 
plored theoretically^ assuming an SU(2) symmetry, as 
appropriate, for example, when the Zeeman energy is 
sufficiently high that only the two-fold pseudospin de- 
generacy remains. In this situation, the graphene FQHE 
problem formally maps into the well studied problem of 
FQHE in GaAs in the zero Zeeman energy limi t 10 ' 11 , with 
the pseudospin of the former playing the role of the spin 
of the latter. In the n = graphene LL, the interaction 
pseudopotential*^ are identical to those in GaAs, so the 
mapping is perfect, and the earlier composite-fermion re- 
sults of Refs. llOllllI carry over to graphene with minimal 
change (spin replaced by pseudospin). In particular, it 
follows that FQHE occurs at fractional fillings given by 



n/(2pn± 1), with pseudospin singlet states at even n and 
pseudospin polarized states at odd n (fully pseudospin 
polarized for n = 1), terminating into a pseudospin sin- 
glet composite fermion Fermi sea in the limit of n — > oo 
(i> = l/2p). The effective interaction in the « / LLs 
in graphene interpolates between those of the |ro|th and 
( | r7. | — l)st LLs for the quadratic dispersion; the composite 
fermion formation^, and therefore the FQHE, are found 
to be almost as strong in the |n| = 1 LL in graphene as 
in the lowest LL of GaAs^. 

It is natural to wonder if FQHE with new structure 
appears in graphene. For this purpose, we explore in this 
work FQHE including the full SU(4) symmetry. By a 
combination of exact diagonalization and the CF theory, 
we find new FQHE states which result from an essential 
interplay between the spin and the pseudospin degrees of 
freedom; such states occur at v = n/(2n + 1) for n > 3. 
For other states, the energy spectrum of the SU(4) prob- 
lem matches with that of the SU(2) problem, although 
the multiplicities are vastly different. We show by ex- 
act diagonalization that the SU(4) symmetry is sponta- 
neously broken at v = | just as at v = 1, i.e., the orbital 
part of the ground state is antisymmetric and the excita- 
tions are skyrmionic. At v — i, | and |, the orbital part 
of ground state is the same as in the SU(2) symmetric 
system, but the state is now a highly degenerate SU(4) 
multiplet. 

The outline of the paper is as follows. We review the 
methods in Sec.|TTJ In Sec. IHII we generalize the compos- 
ite fermion theory to SU(4) electrons, using Fock's con- 
dition and constructing all possible wave functions of in- 
compressible, filled A level states of composite fermions. 
We also indicate which states have analogs in GaAs and 
which do not. Section [IVI is concerned with exact diago- 
nalization results for small systems, which are compared 
to the predictions of the CF theory. Skyrmion-like ex- 
citations are also discussed. In Sec. [V] we obtain ther- 
modynamic energies using the CF theory to determine 
the thermodynamic ground states. Section [VTl discusses 
zero-temperature phase transitions and lists the CF pre- 
dictions for the parameters where the new FQHE states 
should be observable. The paper is concluded in Section 
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In the limit g — > 0, A — > 0, the four vectors 



II. MODEL 

The electronic structure of graphene is well described 
by a tight binding model 3 , which gives a vaulted band 
structure with the valence and conductance bands touch- 
ing in the symmetry protected Dirac (or Fermi) points at 
the six corners of the hexagonal Brillouin zone. The low 
energy states, in particular, occupy double conesi^ with 
apices at the corners. Only two of these six double cones 
are inequivalent, giving rise to a pseudospin degree of 
freedom. Denoting spin by S and pseudospin by P, the 
low energy states are described in the continuum approx- 
imation by an effective Hamiltonian 3 



H = VF{ (a-Ur 



AP z +gn B B-S, (1) 



that acts on a 4-spinor Hilbert space. Here vp s» 10 6 
m/s is the Fermi velocity, II = p+ -A, A is the on-site 
energy difference between the two sublattices, Pi = 1 ® 
!<7, is the pseudospin operator, and S is the spin. Letting 
z = x — iy, measuring distance in units of the magnetic 
length £ B — y/Hc/eB, and using the symmetric gauge 
A = (—By/2,Bx/2), this Hamiltonian is diagonalized 
by the normalized eigenvectors 



(n,m) 
1/2,1/2' 



(n,m) 
1/2,-1/2' 



(n,m) 
-1/2,1/2' 



(n,m) 
-1/2,-1/2 



(8) 



are degenerate, giving rise to an SU(4) internal symme- 
try. 

The SU(4) approximation is valid if Ez = 9HbB and 
A are small in comparison to the scale of the interactions, 
e 2 /(el B ). While A is not easy to estimate, it is generally 
believed to be much smaller than Ez- One can check 
that 



Ez 



d B 



= o.ooie ff V£[rL 



where B[T] is the magnetic field in Tesla. Obviously, 
both e and g depend strongly on the interaction with 
environment in which the graphene sheet is placed. As 
two extremes, one can take e m - m = (e a ; r + esi0 2 )/2 = 2.75 
and e m ax = 15 (the graphite value). On the other hand, 
Zhang et alA find 1.7 < g < 2.0. With these extreme 
values, assuming an extremely high field B = 45 T, we 
get 
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0.032, 
0.208. 
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(g)a s , (5) 



where a 1/2 = , a_ 1/2 = , and r) n<m is the eigen- 

state of the Hamiltonian with quadratic dispersion in LL 
n with angular momentum (L z ) equal to m: 



{-l) n Vn\ 



v / 2 m+1 n(m + n)\ 
The corresponding energies are 



-.z m Ll 



-I-IV4 



(6) 



E nps = sgn(n) 



2hvieB\ 



Ap + gii B Bs. (7) 



Thus the SU(4) symmetric approximation is typically 
valid unless e happens to be close to the high end of 
the estimation. 

The interaction between eigenstates is conveniently de- 
scribed in terms of pseudopotentials^, where the pseu- 
dopotential V m gives the energy of two electrons in 
relative angular momentum m. As a consequence of 
Eqs. ^ and the problem of interacting electrons 
in the nth LL of graphene can be mapped into a prob- 
lem of interacting electrons in rj m states with effective 
pseudopotentials2^ 



y(n)gr. _ 



d 2 k 2w 



F n (k)e- k ~ L m (tf), 



(2tt) 



(9) 



where the form factor F n is 



F (k) 



F n {k) 



L n -i 



, (10) 

The problem of interacting electrons is conveniently 
formulated in the spherical geometry, in which electrons 
move on the surface of a sphere and a radial magnetic 
field is produced by a magnetic monopole of strength 
Q at the centeri 12 ' 15 Here 2Q<f>o is the magnetic flux 
through the surface of the sphere; (fio = hc/e, and 2Q is 
an integer according to Dirac's quantization condition. 16 
For quadratic dispersion the single particle states are 
monopole harmonics 17 YQi m , where I — Q + n is the 
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angular momentum with n = 0, 1, . . . being the LL in- 
dex, m = — I, — I + 1, . . . , I is the z-component of angular 
momentum. The Coulomb interaction is evaluated with 
the chord distance. The analogous solution for carriers 
with linear dispersion, described by the massless Dirac's 
equation, is not known for the spherical geometry. We 
proceed to work in the basis of the monopole harmonics 
YQQm of the lowest LL, while using the effective pseu- 
dopotentials given in Eq. ([9]). Because we have F (k) = 1, 
in the n = Landau level the pseudopotentials derived 
by Fano et alJ^- can be used. In other Landau levels, 
however, we will use the effective pseudopotentials of the 
planar geometry (Eq. ([9])) on the sphere. Such an iden- 
tification is obviously exact in the thermodynamic limit, 
but is usually reasonable also for finite systems. 

We neglect the finite thinkness of the two dimensional 
electron system, which is a much better approximation 
in graphene than in GaAs heterostructures. The neglect 
of Landau level mixing, also made below, is a more tricky 
assumption because both the LL separation (c.f. Eq. ([7|)) 
and the interaction energy scale as V~B, and are of the 
same order. However, their relative strengths can be con- 
trolled by variation of other parameters (Fermi velocity, 
dielectric constant, etc.), and we will assume in the fol- 
lowing that the LL separation is large enough to suppress 
LL mixing. 



III. COMPOSITE FERMION THEORY FOR 
SU(4) ELECTRONS 

We will compare the exact diagonalization results to 
the composite fermion theor y 13 ' 18 , and also use this the- 
ory to explore the thermodynamic limit. (Exact diago- 
nalization can be performed only for very small systems 
for SU(4) symmetry.) The composite fermion (CF) the- 
ory maps the strongly interacting system of electrons in 
a partially filled Landau level to a system of weakly in- 
teracting particles called composite fermions, which are 
bound states of an electron and an even number of quan- 
tized vortices. The mapping consists of attaching 2p 
quantized vortices of the many-body wave function to 
each electron. The quantized vortices produce Berry 
phases that partly cancel the external magnetic field, and 
consequently composite fermions feel a reduced magnetic 
field B* = B — 2pp(f>Q. Composite fermions in magnetic 
field B* fill quantized kinetic energy levels analogous to 
Landau levels, called A levels. If m A levels are filled, 
which corresponds to a filling factor v — 2p m±i > an m ~ 
compressible quantum liquid state results. The FQHE is 
thus explained as the IQHE of composite fermions. 

According to the CF theory, the wave functions for cor- 
related electrons in the lowest Landau level at monopole 
strength Q are given by 

* = Plll$? p $, (11) 



^i = Yl(ujV k ~v jUk ), (12) 

j<k 

where $ is a wave function for N non-interacting elec- 
trons at monopole strength q = Q — p(N — 1), op- 
erator Plll projects a state into the n — LL, and 
u = cos(6»/2)e-^/ 2 and v = sm(6»/2)e^/ 2 . The pro- 
cedure taking us from $ to ^ is called "composite- 
fermionization," which also converts electronic Landau 
levels into A levels of composite fermions. Incompress- 
ible ground states are constructed by completely filling 
the lowest few A levels in one or more spin bands. The 
explicit relation between the monopole strength Q, the 
filling factor v and the particle number N is given below 
in Eq. (fT5j). 

A. Fock's condition 

The above CF construction is valid for the orbital part 
of the wave function independent of whether the Lan- 
dau levels are singly, doubly, or quadruply degenerate, 
i.e., whether the low-energy electronic states have no in- 
ternal (nonorbital) symmetry, have SU(2) symmetry, or 
have SU(4) symmetry, respectively. However, while any 
antisymmetric wave function is valid for a singly degen- 
erate Landau level, only those wave functions are legiti- 
mate for SU(2) which are eigenstates of S 2 and S z . For 
SU(n) electrons, the wave functions are eigenfunctions 
of all the Casimir operators and the generators of the 
Abelian subalgebra of the Lie algebra of SU(n). Among 
such functions the so-called maximal weight states^ sat- 
isfy Fock's cyclic condition 2 ^. The idea will be to con- 
struct valid IQHE wave functions and then show that 
their "composite fermionization" produces FQHE wave 
functions with valid symmetry. 

Fock's cyclic condition tells us which orbital wave 
functions are legitimate for the highest weight state in 
the multiplet [mi, . . . , m n _i] of SU(n). We are using 
the Young tableau notation for the SU(n) multiplets: 
mi > m,2 > . ..m n _i > are integers, rrii being the 
length of the i-th. row of the Young tableau of the repre- 
sentation. (The last to^'s are omitted if zero.) Recall that 
Young tableaux are generated from the direct products of 
the fundamental representation (a single box) with sym- 
metrization in the rows and antisymmetrization in the 
columns. 

Let {a*} be a basis of the (n-dimensional) fundamental 
representation of SU(ti), and let Alt be the number of 
particles assumed to be in the a 1 internal state. Let 

(n max t \ 
n n a ' > ( i3 ) 
t— 1 i— mint / 

where mini = l,maxi = Mi,mhi2 = Mi + l,max2 = 
Mi + M 2 , . . . , and A is the antisymmetrizer. Assume 
^({fj}) is antisymmetric in each subset {min t , . . . , max t } 
of its variables. &'({rj}) is a highest weight state (w.r.t. 
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a choice of positive roots) if and only if it is annihilated f~f~f~f~f~f 

by any attempt to antisymmetrize an electron I of type u f-f-f-f-f-f 

(min u < I < max„) with respect to the electrons of type f-f-f-f-f-f 

t < u > Le - -fHMMMHf- 



ttf 14 4tf 44 



maxj \ 

J2 (M) U({f J }) = 0, (14) " " " (4) 

k— min t / 



$ = $i$2 (15) 



where e is the identity and (k, I) permutes indices k and f f f f f f 

i. 4**m- 

Consequently, any orbital wavefunction f-f-f-f-f-f f-f-f-f-f-f 

ttf 14 ±t 44 

(3,1) 

where $ s 's are Slater determinants such that any state 

(n, to) in <i> s is also filled in (conversely, if (n, m) is 

empty in $ s) then it is also empty in $ s+ i), is a legiti- 

mate highest weight wave function for SU(n). This class f 1 1 f f fl 1 1 f 1 1 f 

includes IQHE states ($s has the lowest l s Landau level f f f f f f f f f f f f 

completely filled; Z s 's are in decreasing order), Fermi sea f4 4 , 1j' 4^ 

states, and single particle, single hole, and particle-hole m 2) 
pair states above the IQHE and Fermi sea states, pro- 
vided that no hole is created in spin state s in a LL that 

is filled in ^s+i, and, conversely, no particle is created in 

spin state s in a LL that is unfilled in f f f f f f 

If state $ satisfies Fock's condition (TJU), so does * f f f f f f f-f-f-f-f-f f-f-f-f-f-f 

given by Eq. (jTT|) . and is therefore a legitimate state. 
Thus composite fermion states with parallel or antipar- 

allel flux attachment can be constructed from electronic (2,1,1) 

Hartree-Fock states as usual. The CF theory for SU(2) 

systems thus generalizes trivially to SU(n) systems. 



ft 14 U 44 



B. New FQHE states 

These observations enable us to construct variational 
wave functions for v = 2p m±i • ^ or m = ^> a sm gl e filled 
A level (in one of the a* spin bands) is the only candidate 
for an incompressible ground state; the SU(4) multiplet is 
[N]. For m = 2, the two A levels are filled with particles 
of identical or different spin state; the multiplet is [N] or 
[N/2, N/2], respectively. As the restriction on the orbital 
part of the CF wave functions is identical to the SU(2) 
case, we expect the latter is prefered in the absence of a 
symmetry-breaking field. 

For m > 3, however, the SU(4) internal degree 
of freedom allows new ground state candidates. Let 
Qii hi hi ^4) denote a trial wave function that fills 
1 1,..., 1 4 A levels of the different single particle spin 
states. As an example, Fig. [TJ depicts all such states for 
to = 4. (For the monopole strength in finite spherical 
models, see Eq. (fl6|) below.) Which of these states will 
be the ground state will be determined only by detailed 
microscopic calculations (to be described below) . A mean 
field approximation, which has been found to be success- 
ful for SU(2) electron o 10 ^ 1 , predicts that the state with 
the lowest "CF kinetic energy" is the ground state. This 
approximation, however, neglects the residual interaction 



H-f-f-f-f -f-f-f-f-f-f WHHHHr HHMMMMf- 

ttf 14 4tf 4^ 

(1,1,1,1) 

FIG. 1: The construction of composite fermion ground states 
for fractions where four A- levels can be filled (P = |, |, etc). 
While the top three states are possible in SU(2) systems, the 
lowest two require an SU(4) symmetric internal degree of free- 
dom. The filled and empty arrows indicate the z-component 
of the spin S z and the pseudospin P z , respectively. This con- 
struction remains valid when S' z and P' z replace S z and P z 
(c.f. Eq. (|19|l ). The top three states also occur in GaAs; the 
bottom three are new. 

between composite fermions, which favors the maximally 
polarized state due to exchange contribution. 

IV. EXACT DIAGONALIZATION RESULTS 

A. Ground states 

Table Q] shows the angular momentum and the SU(4) 
multiplet of the ground state in the |n| = 0, 1 Landau 
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levels at flux values given by 

2Q = 2p(N-l)± ( ^"H, (16) 

where the 1%, h) state is a candidate for the ground 
state at v = 2p "±i with m = ^\ U,. The following ob- 
servations are consistent with the CF theory: 

1. All ground states at the N and Q values satisfying 
Eq. (|16p have rotational symmetry (L = 0), con- 
sistent with the CF theory expectation of incom- 
pressibility for these parameters A and 2Q. These 
systems therefore are finite size representations of 
incompressible states at v — 2p m±i ■ 

2. The filling v = | corresponds to m = 1, a sin- 
gle filled A level (in one of the a* spin bands) is 
the only candidate for an incompressible ground 
state; the SU(4) multiplet is [A], which is consis- 
tent with Table fl] The energy is identical to the 
SU(2) case, but the multiplicity of the ground state 
is different (sa A 3 /6). Analogous to v = 1, the 
ground state is a completely symmetric [A] multi- 
plet, which has associated with it a fully antisym- 
metric orbital wave function, in accordance with 
the Pauli's principle. This is the SU(4) equivalent 
of the well-known ferromagnetic behavior in sys- 
tems with quadratic dispersion relation (GaAs het- 
erostructures and quantum wells) in the vanishing 
Zeeman energy limit. Unlike in the SU(2) sym- 
metric case, however, the spins are not necessarily 
aligned in the same direction; in the [A] multiplet 
spin and pseudospin are equal, but they can take 
all of the possible values^-^: S = P = 0,1, . . . , §■ 
if N is even, and S = P = ±, §, . . . , f if N is odd. 
(The total multiplicity comes from these choices 
combined with the choice of P z and S z .) Any weak 
perturbation that polarizes the spin should also po- 
larize the pseudospin. 

3. For m — 2, two A levels are filled with particles 
of identical or different spin state; the multiplet 
is [A] or [A/2, A/2], respectively. As the orbital 
part of the CF wave functions is identical to the 
SU(2) case, it is no surprise that the latter is pref- 
ered in the absence of a symmetry-breaking field 
(c.f. the f=§,§ data in Table fl}. (The multi- 
plicity of [f , f ] is w A 4 /96.) The ground state 
energy is the same as in the SU(2) symmetric case. 
In Sec. IIIII it was shown that the maximal weight 
state is a product of two Slater determinants just 
as for an SU(2) singlet. Because the SU(2) systems 
already find the lowest energy state in this class, 
the orbital part of the wave function remains the 
same. However, the ground state is neither a spin 
nor a pseudospin singlet; all spin S and pseudospin 
P combinations occur with (S + P) = 0, 2, . . . , y 
for A even, and (S + P) = 1, 3, . . . , f for A odd. 



4. At v = I and v = |, the n = LL ground 
state multiplets are consistent with three and four 
copies of completely filled lowest A levels, respec- 
tively, consistent with the expectation from the 
mean-field model. (Notice that [N/3, N/3, N/3] 
and [A/4, A/4, A/4, A/4] are conjugate represen- 
tations to [A/3] and [0], respectively, which cannot 
be distinguished by their (S, P) quantum numbers.) 
The mean-field expectation is not borne out in the 
I rt I = 1 LL. In Sec. IIIII however, we have seen that 
new composite fermion ground states become possi- 
ble at these fractions; the small system calculations 
we can perform at v = | and v — | are unable to 
determine the true nature of the ground states at 
these fractions. 
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TABLE I: Orbital angular momentum L {n) and S(7(4) multi- 
plet of the ground states of finite systems on a sphere at the 
most prominent fractions in the n = and |n| = 1 Landau 
levels of graphene. The quantity v is the filling factor within 
the Landau level; Q is the monopole strength; N is the num- 
ber of electrons; and D is the dimension of the Hilbert space 
in the L z = S z — Pz = sector. 



B. Charged excitations: SU(4) CF skyrmions 

For SU(2) symmetry and vanishing Zeeman energy, the 
charged excitations at v = 1/3, which map into filling 
factor one of composite fermions, are skyrmions of com- 
posite fermions^. These are analogous to skyrmions of 
electrons- 6 at v = 1. The skyrmions at these fillings 
are obtained as the ground states with the monopole 
strength Q changed by ±1 relative to the Q of the se- 
quence representing the state at v. As Tables [TT] and IIIII 
testify, near v = |, these states belong to an [y, -j] 
or an [-^p-, -^p-] multiplet, depending on whether N is 
even or odd. (The degeneracies of these multiplets are 
w A 4 /96 and w A 4 /48, respectively.) We identify these 
states with SU(4) composite fermion skyrmions. The en- 
ergies of these states are identical to the familiar SU(2) 
skyrmions, as these spinors enforce (cf. Sec. IIII[) exactly 
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the same partition of the wave function into antisym- 
metric factors as the SU(2) singlet (N even) or spin- 
\ (N odd). Consequently, the gaps at v = tttW are 



2p+l 



A { °1 = 0.043(5)e 2 /e/ B and A[ L > = 0.017(3)e7d B 
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TABLE II: Orbital angular momentum L (n) and 517(4) mul- 
tiplet of the ground states of finite systems on the sphere on 
the quasiparticle side of v = 1 and | in the n = and \n\ = 1 
Landau levels of graphene. D is the dimension of the Hilbert 
space in the L z = S z — P z = sector. 



TABLE IV: Coefficients in Eq. {T7J), which reproduces the 
effective interaction for the \n\ = 1 LL of graphene. 



identical to those calculated with SU(2) symmetry£2£: Coefficient 



Value 



co 

Cl 
C2 
C3 
C4 
C5 
C6 



-11.1534 

24.9078 

-18.6461 

6.63657 

-1.221097 

0.112068 

-0.00404269 



regime of validity of the CF theory, we now use the wave 
functions of Eq. (jlljl to study large systems by the Monte 
Carlo method, and obtain thermodynamic limits for the 
energies of various candidate states. 

For the n — Landau level the computation is straight- 
forward. In the \n\ = 1 Landau levels the calculation of 
the energy expectation value by the Monte Carlo method 
requires a knowledge of the real-space interaction corre- 
sponding to the effective pseudopotentials Vm ■ Fol- 
lowing a well-tested procedu re 24 ' 25 , we use a convenient 
form, 



V 


N 


2Q 


D 


L (P) 


5(7(4) m. 


L (1) 


517(4) m. 


1 


4 


4 


117 





[2,2] 





[2,2] 




5 


5 


521 


i 

2 


[3,2] 


i 

2 


[3,2] 




6 


6 


3868 





[3,3] 





[3,3] 




7 


7 


21170 


1 

2 


[4,3] 


1 

2 


[4,3] 




8 


8 


165992 





[4,4] 





[4,4] 


1 

3 


4 


10 


1281 





[2,2] 





[2,2] 




5 


13 


17490 


1 

2 


[3,2] 


1 

2 


[3,2] 




6 


16 


385784 











TABLE III: Orbital angular momentum L (n) and 517(4) mul- 
tiplet of the ground states of finite systems on the sphere on 
the quasihole side of u — 1 and | in the n = and \n\ = 1 
Landau levels of graphene. D is the dimension of the Hilbert 
space in the L z — S z = P z — sector. 

As a technical remark, we note that the Hilbert space 
in the calculations leading to the result of Tables Ul lllll can 
be further reduced by exploiting the existence of a third 
SU(4) generator, traditionally called E 00 , that commutes 
with the Hamiltonian. We identify SU(4) multiplets only 
by their (5, P) quantum number s 21 ' 22 , and do not imple- 
ment this simplification. 



V. THERMODYNAMIC LIMIT 

Because of the four-fold degeneracy of the graphene 
Landau levels, exact diagonalization studies are possi- 
ble only for very small systems. Having established the 



V cS (r) = - + 



M 

£« 

i=0 



(17) 



and fit the coefficients q to reproduce the first M + 1 



pseudopotentials V T 
{Vo 



y(n)cS 



(1)gr ' of Eq. © by 

Y.i <S V eS (zi-z i )\v».. 



{V0,m\V0,: 



(18) 



The coefficients Ci thus obtained are given in Table IIVI 
Our effective interaction produces the Vjj"^ sr ' , . . . , Vq^' 
pseudopotentials exactly; we have checked that the rela- 
tive error in the remaining pseudopotentials is very small 
and does not affect our conclusions. 

We have studied composite fermion wave functions for 
v = 2p m±i ■ Figure [5] shows the energies as a function 
of N for all ground state candidates at these fractions, 
and the thermodynamic limits of the energies are given 
in Table El 

For the sequences v = 2 m+\ > which represent the in- 
tegral quantum Hall effect of 2 CFs (composite fermions 
carrying two vortices), the ground state is consistent with 
the prediction of the mean-field approximation: it is the 
one that best exploits the SU(4) spin bands to minimize 
the CF kinetic energy. The examples are: the state 
(1,1,1) at v = f; the SU(4) singlet state (1,1,1,1) at 
v = |; state (2, 1, 1, 1) at v = and state (2, 2, 1, 1) at 

For the sequence v = 4r ™ +1 , appropriate for 4 CFs, the 
mean-field approximation is not always valid. The fer- 
romagnetic CF ground states are competitive even when 



7 



they do not have the lowest CF kinetic energy, as seen in 
Figure [3] and Table [VI] The ground state at v — ^ is 
found to be fully polarized; the state with the least pos- 
sible CF kinetic energy per particle has slightly higher 
energy. At v = j^, ^- and ^ the states (1,1,1,1), 
(2, 1, 1, 1) and (2, 2, 1, 1), respectively, win by a very small 
margin. Apparently, the gain from CF kinetic energy 
minimization is on the same order as that from exchange 
maximization; which state becomes the ground state is 
determined by their competition. Our results demon- 
strate that the inter-CF interaction are stronger for 4 CFs 
than for 2 CFs. 

For the SU(2) case, it had been found 26 that for v = 
2™, i the model of non-interacting composite fermions 
predicts the ground state quantum numbers correctly, 
which are also confirmed extensively through several ex- 
periments. For v = , on the other hand, a fully spin 
polarized state was found to have the lowest energy in de- 
tailed calculations with the CF theory 26 . The behavior 
is more complicated for the SU(4) symmetry, where, at 
least for some fractions, the fully polarized states are not 
the ground states. 

The data in Tables fVl and I VII include the ground states 
that are available in SU(2) symmetric systems. Our re- 
sults are consistent with those of Park and Jairtii for 
the SU(2) case, although the values differ slightly; the 
current results are more accurate. The extrapolation 
to the theormodynamic limit in our work is based on 
N = 36 — 104 particles, with linear functions fitted on 
the E(l/N) curve under the condition x 2 < N — 2. 

It is in principle possible to apply similar methods 
to the states at v — 2p m-i ' wnere the effective mag- 
netic field B* felt by composite fermions is antiparallel 
to the real external field B. The projection procedure 
for CF wave functions has been elaborated by Moller 
and Simor>2!, but its implementation is impractical for 
the number of A levels and effective monopole strengths 
q that are of interest. (The maximal degree of derivatives 
for negative- B* states is^i 2\q\ +m, while for parallel flux 
attachment it is only to.) Therefore, we do not pursue 
that direction in the present work. 



V 


State 


Energy n — 


Energy n — 1 


tkJ*/N 


3 
7 


(3) 


-0.44178(8) 


-0.4463(14) 


1 




(2,1) 


-0.44704(15) 


-0.4499(7) 


1/3 




(1,1,1) 


-0.45090(3) 


-0.4544(4) 





4 
9 


(4) 


-0.4471(2) 


-0.4527(15) 


3/2 




(3,1) 


-0.4512(1) 


-0.4556(12) 


3/4 




(2,2) 


-0.45244(9) 


-0.4552(6) 


1/2 




(2,1,1) 


-0.45552(12) 


-0.4562(7) 


1/4 




(1,1,1,1) 


-0.45825(5) 


-0.4587(5) 





5 
11 


(5) 


-0.4508(4) 


-0.4545(9) 


2 




(4,1) 


-0.45399(11) 


-0.4555(7) 


6/5 




(3,2) 


-0.45544(16) 


-0.4549(2) 


4/5 




(3,1,1) 


-0.45762(8) 


-0.4560(6) 


3/5 




(2,2,1) 


-0.45888(7) 


-0.4577(8) 


2/5 




(2,1,1,1) 


-0.46084(5) 


-0.4595(9) 


1/5 


6 
13 


(6) 


-0.4527(4) 


-0.4488(19) 


5/2 




(5,1) 


-0.4556(1) 


-0.4531(14) 


5/3 




(4,2) 


-0.4576(4) 


-0.4535(7) 


7/6 




(3,3) 


-0.4579(2) 


-0.4563(6) 


1 




(4,1,1) 


-0.4591(1) 


-0.4571(13) 


1 




(3,2,1) 


-0.4605(1) 


-0.4556(8) 


2/3 




(2,2,2) 


-0.46161(9) 


-0.4585(9) 


1/2 




(3,1,1,1) 


-0.4621(1) 


-0.4595(6) 


1/2 




(2,2,1,1) 


-0.46308(4) 


-0.4604(7) 


1/3 



TABLE V: Energy per particle for incompressible composite 
fermion ground state candidates of different SU(4) symmetries 
a t ^ = 2m+i ' corresponding to the integral quantum Hall 
effect of 2 CFs. For each fraction v the states are listed in 
decreasing order of energy. The last column gives the CF 
kinetic energy (hu)*) per particle in the thermodynamic limit, 
measured relative to the energy of the lowest A level. 



to change to new quantum numbers S' z , P' z by 

's z \ 



u 



p, 



(19) 



VI. QUANTUM PHASE TRANSITIONS 



Either the Zeeman energy E% 
?p 



S, or the 

pscudo-Zeeman energy E% — AP Z , or both, will break 
the SU(4) symmetry. Assuming the symmetry-breaking 
fields are weak, the effect will be to select the most 
favorable member of the ground state multiplet. (No- 
tice that P Zl S z , as well as the third member Eqq of 
the Abelian subalgebra of SU(4) generators, commute 
with the Hamiltonian of Eq. ([1]) in the SU(4) symmetric 
A — > 0,<7 — > limit.) Slightly stronger fields may drive 
zero-temperature phase transitions between the possible 
CF ground states at a fixed filling factor. It is convenient 



with a unitary U to eliminate one kind of Zeeman energy. 
Choosing 



U = 



1 



1 A 



v/TTA 2 \-A 1 



with A = E§ jE% yields 



El' = 



1 



ViTA^ 



-M 



A 



Vl + A 



E P Z ' = 0. 



(20) 



The phase transitions driven by the effective Eg' are 
given in Table IVIII Here we assume a sufficiently weak 
field that does not mix the low-lying states, but only 
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FIG. 2: [Color online] Energy per particle for the novel incompressible composite fermion ground state wave functions for the 
sequence v = m/(2pm + 1). The states are denoted by (li,h, h, U), explained in the text. The upper panels are for the n — 
Landau level, and the lower panels for the \n\ = 1 Landau level. The energies are given in units of e 2 /eZ_g, where Ib = \Jhc/eB 
is the magnetic length. The thermodynamic limits of the the ground state energies are given in Table FVl 
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FIG. 3: [Color online] Same as in Fig. [2] but for the sequence v = m/(4m + 1), corresponding to the integral quantum Hall 
effect of 4 CFs. The thermodynamic limits of the the ground state energies are given in Table IVll 



9 



V 


State 


Energy n — 


Energy n — 1 


hcu*/N 


3 
13 


(2,1) 


-0.34759(4) 


-0.36094(17) 


1/3 




(1,1,1) 


-0.34801(3) 


-0.36112(17) 







(3) 


-0.34822(12) 


-0.36131(29) 


1 


4 
17 


(2,2) 


-0.35040(8) 


-0.3642(2) 


1/2 




(3,1) 


-0.35062(7) 


-0.3648(2) 


3/4 




(2,1,1) 


-0.35081(3) 


-0.3643(2) 


1/4 




(4) 


-0.35106(5) 


-0.3646(3) 


3/2 




(1,1,1,1) 


-0.35114(3) 


-0.3650(1) 





5 

1>1 


(3,2) 


-0.35212(5) 


-0.3652(3) 


4/5 




(4,1) 


-0.35238(9) 


-0.3658(4) 


6/5 




(3,1,1) 


-0.35257(5) 


-0.3659(3) 


3/5 




(2,2,1) 


-0.35257(5) 


-0.3663(1) 


2/5 




(5) 


-0.35266(10) 


-0.3677(3) 


2 




(2,1,1,1) 


-0.35283(6) 


-0.3666(3) 


1/5 


6 
26 


(3,3) 


-0.35328(2) 


-0.3661(2) 


1 




(4,2) 


-0.35337(6) 


-0.3666(2) 


7/6 




(5,1) 


-0.35354(9) 


-0.3668(4) 


5/3 




(4,1,1) 


-0.35366(5) 


-0.3665(1) 


1 




(3,2,1) 


-0.35375(3) 


-0.3667(1) 


2/3 




(2,2,2) 


-0.35387(7) 


-0.3670(2) 


1/2 




(3,1,1,1) 


-0.35393(3) 


-0.3669(1) 


1/2 




(6) 


-0.35394(18) 


-0.3663(6) 


5/2 




(2,2,1,1) 


-0.35405(4) 


-0.3674(2) 


1/3 



TABLE VI: Same as in Table |V] for v = correspond- 
ing to the integral quantum Hall effect of 4 CFs. Energies per 
particle for incompressible composite fermion ground state 
candidates with different SU(4) symmetries are shown in de- 
creasing order. For each fraction v the states are in decreasing 
order of energy. 

selects the most favorable member of the ground state 
multiplet. In state (Zi, 1%, I3, Z4) the with li > h+i, this 
means filling h + h A levels of the favorable S' spin, and 
l 3 + I4 levels of the unfavorable S' spin. This results in 
an effective Zeeman energy difference per particle. No- 
tice that no value of E% and E§ will drive the system 
at v — m/(2pm + 1) to a completely antisymmetric or- 
bital state (m). This is a consequence of our freedom to 
choose a basis P' z , S' z in the Cartan subalgebra of SU(4) 
(Eq. (|19p ). which eliminates one kind of Zeeman energy; 
the states where the two unfavorable bands of the S' spin 
are emptied will already have the lowest possible effective 
Zeeman energy E% . 

Tilting the magnetic field is often used as a means 
to tunc the effective Zeeman energy. Unlike in 
GaAs/AlGaAs heterostructures, the in- plane megnetic 
field is unlikely to change the single-particle states in 
any significant manner, because the transverse thick- 
ness of the two dimensional electron system is graphene 



(~ 2A) is much smaller than the typical magnetic lengths 
(^100A). The most favorable parameter space for the ob- 
servation of these transitions occurs for very low values 
of the dielectric constant e and the sublattice asymmetry 
A. The value of e depends on the interaction with the 
substrate on which the graphene sheet is placed. 



V 


Transition 


( E * \ 


Change of (S') 


3 
7 


(1,1,1) -(2,1) 


0.0116(6) 


11 

6 2 


4 
9 


(1,1, 1,1) -(2, 1,1) 
(2, 1,1) -(2, 2) 


0.0109(7) 
0.0123(8) 


11 

4 2 


5 
11 


(2, 1,1,1) -(2, 2,1) 
(2, 2,1) -(3, 2) 


0.0098(6) 
0.0172(8) 


13 
10 10 

3 1 
1(1 2 


6 
13 


(2,2,1,1)^(3,2,1) 
(3, 2,1) -(3, 3) 


0.0076(5) 
0.016(2) 


11 

6 3 
1 1 
3 2 



TABLE VII: Zero temperature phase transitions driven by the 
effective Zeeman energy in the n = Landau level. The last 
column gives the expectation value of the transformed spin 
S' (Eq. JHJ), which coincides with the ordinary spin in the 
case of perfect sublattice symmetry (A = 0). The quantity 

(jSg 1 1 (e 2 /eZs)^ give the value of the parameter where the 

transition is predicted to take place. 



VII. CONCLUSION 

We have used a combination of exact diagonalization 
and the CF theory to identify a large range of possi- 
ble FQHE states in graphene allowed by SU(4) symme- 
try, and shown that new states, which have no analog in 
GaAs, can occur at filling factors v — 2p " ±1 for m > 3 
in the \n\ = and \n\ = 1 Landau levels of graphene. 
For 2 CFs, the ground states are those for which the com- 
posite fermion kinetic energy is minimum; these states 
spread out maximally in SU(4) spin space. For 4 CFs 
the fully polarized ground state wins, with the excep- 
tion of v = A, where very compact SU(4) singlet struc- 
tures are possible. We have also estimated parameter 
regimes where these states should occur; zero tempera- 
ture phase transitions can be driven by variation of an 
effective Zeeman energy that accounts for the combined 
effect of the Zeeman energy and the sublattice symmetry 
breaking field. At v — | we found SU(4) skyrmions to 
be the lowest energy charged excitations. 
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